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ADAPTIVE  ANTENNA  SUBARRAYING 
USING  A  WEIGHTED  BUTLER  MATRIX 


INTRODUCTION 

A  fully  adaptive  array  is  one  in  which  every  element  of  the  array  is  individually  controlled  adap¬ 
tively.  In  theory,  this  provides  the  necessary  degrees  of  freedom  to  lower  all  of  the  deterministic 
sidelobes  to  any  arbitrary  level.  A  partially  adaptive  array  is  one  in  which  elements  are  controlled  in 
groups  (the  subarray  approach)  or  in  which  only  certain  elements  called  auxiliary  elements  are  made 
controllable. 


The  partially  adaptive  array  has  not  been  extensively  studied  as  evidenced  by  the  lack  of  refer¬ 
ences  available  in  the  literature  [1-5].  Chapman  [3]  initiated  one  of  the  earliest  studies  (circa  1974) 
into  exploring  how  to  effectively  combine  the  N  elements  of  an  entire  planar  antenna  array  into  a  col¬ 
lection  of  L  subarrays.  He  developed  the  row-column  precision  antenna  array  (RCPAA)  in  which  each 
element  signal  is  split  into  two  paths:  a  row  path  and  a  column  path.  All  the  elements  of  a  given  row  or 
column  are  added  together,  and  all  the  row  outputs  and  column  outputs  are  then  adaptively  combined. 
Morgan  [2]  studied  the  technique  of  just  using  auxiliary  elements  to  perform  the  adaptation.  He 
showed  that  the  correllation  coefficient  between  the  adaptive  element  spatial  vectors  plays  a  key  role  in 
characterizing  the  performance  of  a  partially  adaptive  array.  A  significant  difference  in  array  perfor¬ 
mance  can  occur  if  the  array  elements  selected  for  adaptive  control  are  not  properly  chosen. 

Obviously,  the  fully  adaptive  configuration  is  preferred  since  it  offers  the  most  control  over  the 
response  of  the  array.  However,  the  typical  array  may  have  many  elements.  This  poses  several 
immediate  problems.  Implementing  an  adaptive  array  with  many  degrees  of  freedom  can  have  consid¬ 
erable  impact  on  the  total  cost  of  the  array.  This  alone,  however,  is  not  sufficient  grounds  for  rejecting 
the  idea.  Processor  implementation  does,  however,  pose  a  more  serious  problem  for  a  system  of  this 
size.  In  order  to  appreciate  this,  it  is  necessary  to  elaborate  somewhat  upon  the  ways  in  which  such  a 
processor  can  be  implemented. 


For  an  adaptive  processor  with  N  degrees  of  freedom,  each  of  the  N  input  channels  is  multiplied 
by  a  complex  weight  and  summed  to  form  the  output.  The  adaptive  weights  are  determined  according 
to  a  control  law  [6] 


MW  -  MS, 


(1) 


J 

where  M  is  an  N  x  N  matrix  of  the  cross-covariances  of  the  signals  in  the  N  channels,  W  is  an  N- 
element  vector  of  the  weights,  n  is  an  arbitrary  nonzero  constant,  and  S  is  a  vector  representation  of 
the  desired  signal  in  each  channel.  Solving  for  the  weights  is  equivalent  to  solving  a  set  of  N  simul¬ 
taneous  linear  equations.  The  most  straightforward  way  in  which  to  implement  a  processor  to  do  this  is 
to  use  a  form  of  the  Howells  control  loop  in  analog  hardware  [7],  or  the  Widrow  LMS  algorithm  (8],  or 
one  of  its  derivatives  in  a  digital  machine.  Both  are  realizations  of  the  "steepest  descent"  method,  an 
iterative  technique  employed  in  the  minimization  of  a  function.  This  is  a  reliable  method  of  solution  in 
this  application  since  it  can  be  shown  that  the  surface  described  by  the  output  residue  as  a  function  of 
IF  is  a  concave  quadratic  hypersurface  and,  as  such,  contains  only  one  minimum.  The  problem  with 
this  approach  is  that  the  convergence  rate  depends  upon  the  particular  physical  configuration  of  the 
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interrerence  sources  and  the  array.  Mathematically,  the  transient  response  oT  the  processor  is  a  func¬ 
tion  of  the  eigenvalues  of  the  covariance  matrix  M.  The  spread  in  the  amplitude  of  the  eigenvalue  can 
be  considerable  and,  since  the  particular  configuration  of  noise  sources  to  be  encountered  is  never 
known  a  priori,  the  resultant  variance  in  settling  time  can  be  large,  particularly  if  the  order  of  A#  is 
large,  i.e.,  a  targe  number  of  degrees  of  freedom  are  used. 

Reed  et  al.  [9]  have  proposed  to  solve  directly  for  the  weights  by  inverting  an  estimate  of  the  true 
covariance  matrix,  thereby  eliminating  the  uncertainty  of  the  convergence  rate  of  the  iterative  method. 
This  method,  too,  is  subject  to  limits  related  to  the  order  of  the  matrix,  M.  The  number  of  arithmetic 
steps  required  to  obtain  the  solution  by  this  method  is  proportional  to  the  cube  of  the  order,  N.  This, 
coupled  with  the  requirement  for  high  numerical  accuracy,  puts  an  upper  bound  on  N  which  is  related 
to  the  computation  time  available  and  on  the  cost  that  can  be  tolerated. 

In  addition,  if  only  auxiliary  elements  are  chosen  as  inputs  to  the  sidelobe  canceller,  jammer  nul¬ 
ling  is  effected  by  the  bandwidth-aperture  product;  i.e.,  if  the  number  of  antenna  elements  is  held  con¬ 
stant,  as  the  bandwidth  of  the  jammer  or  the  size  of  the  antenna  aperture  in  terms  of  wavelengths  or 
the  spacing  of  the  auxiliary  elements  is  made  larger,  the  nulling  of  the  jammer  is  degraded.  This  degra¬ 
dation  results  because  as  auxiliary  antenna  elements  are  spaced  farther  apart,  the  input  noise  signal  with 
a  spread  spectrum  becomes  more  uncorrelated  between  elements.  However,  it  is  the  correlation 
between  the  signals  on  the  respective  antenna  elements  that  allows  the  adaptive  array  to  sense  the 
direction  of  arrival  of  an  undesirable  signal  and  then  place  a  null  in  the  antenna  pattern  in  the  direction 
of  that  undesirable  signal. 


All  of  these  considerations  point  toward  the  desirability  of  reducing  the  dimensionality  of  the  pro¬ 
cessor  while  maintaining  as  much  control  as  possible  over  the  size  of  a  given  aperture.  Reduction  of 
the  dimensionality  of  the  processor  is  achieved  through  subarraying  as  illustrated  in  Fig.  1.  Here,  the 
dimensionality  is  reduced  from  N  to  £,  where  L<N ,  by  using  a  transformation  matrix  H  which  is  a 
L  x  N  matrix.  It  is  not  a  trivial  problem  to  specify  an  effective  H  matrix.  If  the  H  matrix  is  specified 
improperly,  the  overall  signal-to-noise  ratio  can  degrade  significantly  as  a  function  of  the  direction  of 
arrival  angles  of  the  jamming. 
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Adaptive  subarray  configuration 


This  report  describes  an  effective  procedure  for  specifying  the  H  matrix.  However,  the  H  matrix 
will  not  be  a  constant  as  in  the  RCPAA  (3]  approach  but  it  will  be  a  function  of  the  noise  environment. 
In  a  sense,  the  H  matrix  will  adapt.  We  will  see  in  the  next  section  that  an  effective  choice  for  H  is  a 
partition  of  a  weighted  Butler  matrix. 

The  next  section  describes  the  algorithm  of  specifying  H  in  detail.  Then  results  are  presented  to 
support  utility  of  the  subarray  algorithm.  Specifically  it  is  shown  that  the  maximum  signal-to-noise  ratio 
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possible  using  no  subarraying  suffers  a  small  degradation  if  an  adaptive  weighted  Butler  matrix  subarray- 
ing  scheme  is  used.  Next,  the  advantages  of  using  this  subarraying  scheme  are  presented. 

it  should  be  noted  that  in  this  treatment  of  subarraying  for  simplicity's  sake  we  deal  only  with 
linear  antenna  arrays.  Later  reports  will  consider  planar  arrays. 


THE  WEIGHTED  BUTLER  MATRIX  TRANSFORMATION 


Let  the  transformation,  //,  seen  in  Fig.  I  be  a  Butler  matrix.  A  Butler  matrix,  B,  has  the  form 
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where 


rN 


-j— 
-  1  N 


(3) 


and  j  -  >/— T.  Note  that  in  this  case  L  —  N.  If  v  signifies  the  input  vector  and  u  the  output  vector 

then  u  —  B v  where  u  —  («|,  .. .  ,  u#)T  and  v  «■  (vj . vN)r.  (Note  that  we  use  u  as  the  output  of 

the  transformation  instead  of  xN  «•  (xt,  ....  xL)T as  depicted  in  Fig.  1.  We  do  this  because  we  will  use 

a  subset  of  the  uk,  k  -  1,  2 . N  outputs  as  inputs  to  the  adaptive  processor.  These  inputs  will  be 

denoted  by  5.)  An  output,  uk,  of  the  transformation  has  the  form 


(4) 


The  form  of  Eq.  (4)  is  recognized  as  that  of  a  discrete  Fourier  transform.  It  transforms  v„  which 
is  a  function  of  the  position  on  the  array  (array  space)  indicated  by  n ,  into  angular  beam  space.  Thus, 
if  there  is  a  jammer  located  0  degrees  off-boresight  such  that 

v-  (1,  eJ*.  e1* . e(N-uj*)T  (5) 


where  —  (2ird/k)  sin®,  d  is  the  element  separation,  and  A  is  the  wavelength,  then  the  transformation 
indicated  by  Eq.  (4)  places  most  of  the  jammer  power  into  the  transformation  output,  uk ,  where  k  is 
the  integer  value  that  makes  \(2vk/N)  —  dt  a  minimum. 


To  illustrate  this  further,  consider  the  output  response  of  uk  vs  4.  We  can  show  using  Eqs.  (4) 
and  (S)  that 


,in>*^ 


;  a  - 


2wk 
N  ' 


(6) 


The  output  %  as  a  function  of  d  is  shown  in  Fig.  2.  Observe  that  uk  is  relatively  large  if  4  is  dose  to 
(2 wk/N)  but  decreases  rapidly  as  1(2 wk/N)  -  dl  increases.  Hence,  the  transformation  by  a  Butler 
matrix  tends  to  orthogonalize  the  outputs  of  the  matrix  transformation,  i.e., 
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(7) 


for  /  k.  This  approximation  tends  to  become  better  as  1/  —  k\  increases. 

As  has  been  shown  by  past  researchers  [9- IS],  the  orthogonalization  (For  example  Gram-Schmidt 
orthogonalization  developed  by  Lewis  and  Kretschmer)  of  the  array  inputs  by  preprocessing  enhances 
the  convergence  rate  of  obtaining  the  optimal  adaptive  weights.  The  convergence  rate  using  preprocess¬ 
ing  is  essentially  independent  of  the  eigenvalues  of  the  input  covariance  matrix  and  hence,  does  not 
suffer  the  limitation  of  the  control  loop  implementations.  Hence,  the  settling  times  of  the  weights  are 
almost  as  fast  as  those  obtained  by  using  direct  matrix  inversion  [9]. 

We  can  extend  the  concept  of  using  a  Butler  matrix  transformation  to  that  of  using  a  weighted 
Butler  matrix  where 
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Note  that  we  can  write  H  as 


H-  BA 


(9) 


where  A  is  a  diagonal  matrix  with  diagonal  elements  ak,  k  —  1,  ,  N.  The  reason  for  using  the 

weights,  ak,  is  to  control  the  near-in  sidelobe  levels  of  the  response  of  uk  vs  However,  by  lowering 
the  near-in  sidelobe  levels  we  also  enlarge  the  beamwidth  of  the  response  of  uk,  k  —  1,  ....  N  vs 
If  this  occurs,  it  may  then  be  desirable  to  beamform  as  an  output  every  other  uk,  or  every  third  u*,  etc. 
Hence,  by  controlling  these  sidelobes  we  can  make  adjacent  outputs  more  uncorrelated  and  therefore 
enhance  the  convergence  rate.  However,  by  doing  this  we  may  lose  some  desired  signal  power  ratio 
because  now  H  may  not  be  unitary  transformation. 


4 


NRL  REPORT  8636 


For  the  subarraying  technique  to  be  presented,  we  only  use  L  of  the  Butler  transformed  outputs 
where  L  <  N.  We  choose  the  L  outputs  as  follows  (see  Fig.  3):  Let  (xk),  k  -  1 . L  be  .the  sub¬ 

set  of  {«/*)  k  —  1,  ....  N  that  is  chosen  as  inputs  to  the  adaptive  processor  such  that  X|  is  always 
chosen  as  one  of  the  L  outputs.  It  can  be  shown  directly  that 

N 

*i  -  £  ak\k.  (10) 

Thus  we  use  X|  as  the  main  antenna  channel  where  the  weights  ak,  k  —  1 . N  can  also  be  specified 

to  control  the  steering,  sidelobe  level,  beamwidth,  ect.  of  the  antenna.  Therefore,  we  see  that  the 
adaptive  antenna  is  configured  as  a  sidelobe  canceller.  After  multiplication  by  the  weighted  Butler 

matrix  the,  uk ,  k  —  2 . /V  are  ordered  from  largest  to  smallest  based  on  their  individual  powers. 

Next  the  L  -  1,  «*,  k  »  2,  ... ,  L  voltages  with  the  largest  powers  are  selected  as  the  inputs,  xk, 
k  »»  2 . L  to  the  sidelobe  canceller. 
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Fig.  3  —  Functional  block  diagram  of 
subarraying  algorithm 


The  largest  powers  of  the  outputs  were  chosen  based  on  the  analysis  rendered  in  Appendixes  A 
and  B  and  the  next  section.  There  we  demonstrate  empirically  (and  analytically  in  a  specific  case)  that 
the  sidelobe  canceller  suffers  minimum  loss  in  signal-to-noise  ratio  if  the  largest  powers  of  the  outputs 
of  the  weighted  Butler  matrix  transformation  are  selected. 


It  is  evident  from  the  procedure  defined  that  the  number,  L,  of  subarray  outputs  can  be  a  vari¬ 
able.  For  example,  if  there  are  only  Nj  spatially  distinct  narrowband  jammers  then  L  —  Nj.  Also,  if 

we  set  a  power  threshold  such  that  only  those  outputs,  uk,  k  -  2 . Nthat  exceed  this  threshold  are 

chosen,  then  L  will  be  a  function  of  the  threshold,  the  noise  environment,  and  the  sidelobe  levels  that 
result  from  weighting  the  Butler  matrix. 


RESULTS 

One  of  the  costs  of  subarraying  is  the  maximum  signal-to-noise  ratio  that  is  possible,  (S/N)wh, 
degrades  from  the  maximum  signal-to-noise  ratio,  (S/N),  that  is  possible  with  no  subarraying.  We 
define  this  loss  to  be 
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<^«.b)dS-  Nlog 


(sm* 


where  I  ■  I  denotes  absolute  value.  We  demonstrate  that  using  the  subarraying  scheme  described  in  the 
previous  section  for  the  cases  considered  results  in  a  small  loss. 

Let  us  first  consider  the  following  optimization  problem:  given  that  there  is  a  narrowband  jammer 
9j  degrees  off-boresight  with  jamming  power  &}  normalized  to  the  internal  noise  power,  what  2  x  N 
subarray  transformation  using  a  partition  of  the  Butler  matrix  with  uniform  weights  results  in  the 
minimum  subarray  (S/A/)  loss,  Lsub? 

In  more  mathematical  terms,  if 


where 


H-  B„ 

1  (1  1  1 


V77  1  e~j*k  e~V*k  '  •  ■  e-w-M**r 


i  2ir  . 


1.2 . N-  1. 


then  what  value  of  the  integer  k  results  in  the  minimum  loss  as  defined  by  Eq.  (ID?  Note  B„  is  a  par¬ 
tition  of  the  Butler  matrix,  B. 

For  this  problem,  expressions  for  (S/A/),„b  and  (S/AO  are  derived  in  Appendix  A  (also  in  this 
appendix  an  equivalent  transform  to  the  one  seen  in  Fig.  1  is  derived  and  it  is  shown  in  general  that  for 
any  subarray  transformation,  (S/N)^  can  be  found  directly  from  the  formulated  expression  of 
(S/AO).  These  expressions  were  used  in  Appendix  B  to  show  that  k  should  be  chosen  to  minimize 
I*  -  **l,  where  *  -  (2 nd/\)  sing.  Hence,  should  be  selected  as  close  to  the  value  of  *  as  possi¬ 
ble  to  obtain  the  minimum  S/ A/ loss. 

It  is  also  shown  in  Appendix  B  that  if  a  2  x  N  weighted  Butler  matrix  is  used  such  that  |*  —  **| 
is  a  minimum,  and  No]  »  1  then  the  maximum  S/ N  loss  in  dB  is  approximately 


where 


(Lj3bAX)dB-  10.7161  |f(*)  |  |J  dB, 

f<*f-  is  •**■'** 


and  It'll  denotes  complex  magnitude.  Note  that  the  quiescent  (no  jamming)  antenna  pattern  can  be 
defined  in  terms  of  a*,  k  —  1,  ....  N and  is  equal  to  Ilf  (a)|(2  if  we  normalize  all  antenna  gains  to  the 
mainbeam  gain  of  a  uniformly  weighted  array,  a  -  (2 wd/k)  sin/8  and  0  is  an  arbitrary  angle  off- 
boresight.  We  see  from  Eq.  (IS)  that  the  iS/N)  loss  in  d&  will  be  small  if  the  jammer  is  located  in  the 
sidelobes  and  is  proportional  to  the  normalized  antenna  gain  in  the  direction  of  the  jammer.  For  exam¬ 
ple,  if  the  quiescent  normalized  antenna  gain  in  the  direction  of  the  jammer  is  -20  dB,  then  the  iS/N) 
loss  calculated  using  Eq.  (IS)  is  approximately  no  greater  than  0.1  dB. 

To  further  illustrate  that  choosing  the  k  that  minimizes  I*  -  **|  results  in  the  minimum  SIN  km 
we  have  plotted  this  loss  vs  k  and  ***-1,2,  ....  N~  1  fnr  the  following  cases: 

1.  A/-  »,  L-l,  gy-  10°  (Fig.  4), 

2.  N  -  16.  L  -  2,  9j  -  10*  (Fig.  5), 

3.  A-  16.  L  -  2.  9j  -45*  (Fig.  6). 
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In  these  cases,  a  uniformly  weighted  Butler  matrix  is  used.  We  see  from  Figs.  4,  S,  and  6  that  indeed 
the  minimum  (S/N)  loss  occurs  when  \</>  -  >lik I  is  a  minimum.  For  cases  1  and  2,  <f>  —  31  “and  for 
case  3,  <j>  -  127°.  Note  from  comparing  the  curves  of  Fig.  5  with  Fig.  6,  that  the  minimum  loss  is 
greater  as  the  jammer’s  direction  of  arrival  (DOA)  moves  closer  to  boresight.  This  observation  is 
predicted  by  Eq.  (15). 


Now  consider  the  case  where  there  are  two  jammers  and  the  number  of  antenna  elements  is  36. 
Both  jammers  have  <rj  *20  dB  power  and  the  first  jammer  is  0i  —  25  °  off-boresight.  Let  us  vary  the 
DOA,  02i  of  the  second  jammer  from  0°  to  90°  off-boresight  and  subarray  using  a  3  x  Af  uniformly 
weighted  partitioned  Butler  matrix,  such  that 


B, 


1  1  1 

~V*kk 

1  e  e 


i  ! 

-w-u/aJ 

e 

e 


(17) 


In  the  above  equation  (2it/N)k\,  *lik}  -  (lir/N)k2  and  k\,  k2  are  chosen  to  minimize 

l<£i  -  0^1  and  \<f>2  -  1  (<^j  -  31  <f>2  —  it  sin02)  such  that  k{  &  k2  (note  if  —  k2  then  the  H 

transformation  is  of  rank  2  which  renders  the  calculation  of  the  ( S/N )  loss  impossible  because  of  the 
necessity  of  inverting  a  singular  matrix;  see  Appendix  A,  Eq.  (A5)).  If  it  results  that  k\  -  k2,  then  we 
choose  k2 *»  k\  +  1  or  k2  -  -  I  such  that  k2  ^  0  MOD  (/V). 

The  resultant  (S/S)  loss  is  graphed  in  Fig.  7.  We  have  not  attempted  to  plot  L, 1B  for  every  92 
but  have  plotted  L, dB  in  increments  of  A92  -  l°-  First  note  that  in  the  sidelobes  for  this  case  that  the 
average  loss  is  small  (=0.008  dB).  Secondly,  the  loss  becomes  larger  as  the  second  jammer  moves 
through  the  boresight  of  the  antenna.  This  is  not  unexpected  because  the  null  in  the  direction  of  the 
jammer  also  nulls  the  desired  signal  as  the  null  is  placed  closer  to  the  main  beam.  Hence  a  loss  of 
(S/N)  performance  occurs. 

The  undulating  nature  of  the  subarray  loss  is  due  to  the  undulating  nature  of  l4>2  -  I.  A  plot 
of  1^2  ~  appears  in  Fig.  8.  As  l$2  -  2 1  increases,  the  subarray  loss  becomes  larger  because  the 

phase,  $2*  is  n°l  centered  on  one  of  the  discrete  angles,  (2ir/AD  •  k,  k  -  1,  2 . N  —  I.  Also  note 
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in  Fig.  7  that  as  the  second  jammer  angle  passes  thru  the  first  jammer  angle  (02  ~  25*7,  there  is  at  first 
a  very  small  (S/N)  loss  (~  0.0016  dB)  but  then  a  noticeably  large  (S/N)  loss  (==  0.19  dB).  We  will 
explain  this  phenomenon  in  a  future  report. 

ADVANTAGES  AND  FUTURE  STUDIES 

Some  of  the  major  advantages  of  subarraying  as  described  in  the  preceding  sections  are  that: 

1.  it  reduces  the  dimensionality  of  the  sidelobe  canceller  processor  which  reduces  the 
amount  of  necessary  hardware, 

2.  it  effectively  cancels  sidelobe  jamming  without  significant  degradation  of  the  output 
signal-to-noise  ratio. 
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3.  reduction  of  the  dimensionality  of  the  sidelobe  canceller  decreases  the  convergence  time 
of  the  sidelobe  canceller  (for  the  matrix  inversion  algorithm  [1]  convergence  rate  is  pro¬ 
portional  to  the  dimensionality), 

4.  convergence  time  of  the  sidelobe  canceller  is  independent  of  the  size  of  the  array;  it  is 
proportional  to  the  number  of  jammers, 

5.  orthogonalization  of  the  channels  allows  the  theoretical  convergence  time  to  be  achieved 
since  the  loops  are  now  almost  independent  and  do  not  interact, 

6.  the  effect  of  the  bandwidth-aperture  product  is  reduced  since  all  the  elements  of  the 
antenna  are  used  to  form  the  reduced  number  of  inputs  to  the  sidelobe  canceller. 

The  third  advantage  mentioned  above  results  because  if  the  direct  matrix  inversion  algorithm  (91  is 
used  to  compute  the  optimal  weights,  then  the  average  number  of  samples  of  each  subarray  output 
required,  such  that  the  average  output  signal-to-noise  ratio  is  within  3  dB  of  the  optimum  subarrayed 
signal-to-noise  ratio,  is  equal  to  1L  where  L  is  the  number  of  subarrayed  outputs.  If  we  do  not  use 
subarraying  but  use  the  entire  array  then  the  number  of  samples  of  each  array  output  needed  is  IN. 
Obviously  if  L  «  N,  then  the  number  of  samples  needed  when  subarraying  is  much  less  when  not 
subarraying.  Hence,  the  convergence  rate  will  be  faster.  (Note  that  it  is  as  assumed  that  the  samples 
are  from  time  independent  zero-mean  Gaussian  processes.) 

The  fourth  advantage  is  true  under  the  assumption  that  we  are  thresholding  the  powers  of  the 
individual  output  channels  of  the  weighted  Butler  transformation. 

Future  studies  will  be  directed  into  investigating  the  effect  of  random  errors  in  the  antenna  ele¬ 
ments.  We  have  assumed  in  our  analysis  that  all  of  the  antenna  elements  are  identical.  In  practical 
cases  this  is  not  true.  There  will  be  small  random  phasor  errors  which  multiply  the  outputs  of  the  indi¬ 
vidual  antenna  elements.  These  random  errors  will  degrade  S/N  performance  and  upper  bound  the 
maximum  S/N. 

A  second  study  will  be  undertaken  to  investigate  the  effects  of  the  jammer  bandwidth-antenna 
aperture  product  on  the  subarry's  performance.  Since  this  product  inherently  limits  the  performance 
of  a  totally  adaptive  array,  it  may  or  may  not  further  limit  the  performance  of  a  subarrayed  adaptive 
antenna  where  we  subarray  using  a  weighted  Butler  matrix. 
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Appendix  A 

EQUIVALENT  TRANSFORMATIONS  FOR  SUBARRAYS 

In  this  appendix,  we  prove  that  the  subarraying  transformation  as  seen  in  Fig.  A1  is  equivalent  to 
the  transformation  seen  in  Fig.  A2  in  that  the  resultant  maximum  signal-to-noise  ratios  are  the  same. 
The  //transformation  seen  in  Fig.  A1  is  an  L  x  N matrix  and  has  rank  L. 

The  transformation  viewed  in  Fig.  A 1  works  as  follows:  the  internal  (thermal  noise)  and  external 
inputs  add  at  the  front  end  (the  antenna  elements)  after  which  an  L  x  N  matrix  transformation  reduces 
the  N  inputs  to  L  outputs  (L  <  N).  Finally,  these  L  outputs  are  optimally  weighted  such  that  the 
maximum  signal-to-noise  ratio  is  achieved. 

The  transformation  as  seen  in  Fig.  A2  is  not  actually  realizable  because  here  we  are  separating  the 
internal  noises  from  the  external  noises  and  operating  on  them  separately.  This  transformation  works 
as  follows:  the  N  external  inputs  are  multiplied  by  an  N  x  N  matrix  P  where  P  —  H'{HHr)~x  H  and  t 
denotes  conjugate  transpose.  The  outputs  of  this  transformation  are  then  added  to  the  internal  inputs 
(thermal  noise  of  the  same  rms  power  as  the  thermal  noise  seen  in  Fig.  Al).  Finally,  the  resultant  N 
outputs  are  optimally  weighted. 


Let  X  represent  the  noise  vector  of  N  external  antenna  inputs,  rj  represent  the  vector  of  N  inter¬ 
nal  noise  sources,  and  V  represent  the  total  noise  vector.  Obviously 

Y-X  +  ij,  (Al) 

and  _  _ 

E\  YW  -  Em  +  Efntjl,  (A2) 

assuming  the  internal  noise  processes  are  zero  mean  and  independent  of  the  external  noises.  Further¬ 
more,  let  the  internal  noise  sources  be  statistically  independent  and  identically  distributed  with  noise 
power  equal  to  <r$.  Thus 
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i| 

>  i 
if 


£,{nn,l  —  «rii 

(A3) 

where  /  is  the  /V  x  N  identity  matrix. 

If  we  set  M-  EfPPl  and  Mj 

-  £{ JFjPl,  then 

M-*$l  +  Mj. 

(A4) 

If  5  represents  the  input  signal  vector,  then  it  can  be  shown  [A  1_1  that  the  maximum  signal-to-noise 
ratio,  (S/N)\  that  is  possible  by  properly  weighting  the  outputs,  7,  of  the  subarray  transformation  as 
seen  in  Fig.  A1  is 

S 

N 

-  S'H'i  HMHt)  ~ 1 HS. 

i 

(A5) 

We  can  also  show  that  the  maximum  signal-to-noise  ratio,  ( S/N)2 ,  that  is  possible  by  properly  weight¬ 
ing  the  output,  Z,  as  seen  in  Fig.  A2  is 

S 

N 

l 

(A6) 

where 

Sp  -  PS. 

(A7) 

Mr 

-  <7q/  +  PMjP1. 

(A8) 

and 

P- 

(A9) 

We  will  show  that  (S/AO*  *»  (S/W)2. 

To  this  end,  let  us  normalize  the  covariance  matrices  to  the  internal  noise  power  such  that 

and 


It  is  easy  to  show  that 


and 


£ 

N 


where 

Mr 

We  write  ( S/N)2  in  a  more  expanded  form 

S 


lit  -  I  +  Mj. 

(A  10) 

Mj-  "T  Mj. 

(All) 

~  S'H'iHMHT'HS. 

*6 

(A12) 

it  si>Mp'sP’ 

(A13) 

/  +  PMjP'. 

(A  14) 

-  S'P'M-'PS. 

t  <r( 

(A15) 
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We  prove  that  ( S/N)\  —  ( S/S)2  by  showing  that 

-  P'M~'P.  (A16) 

Before  we  do  this,  let  us  state  that  it  can  be  shown  that  P  is  an  idempotent  matrix,  i.e.,  Z*2  -  P  and 
therefore  singular  for  /V  >  2.  In  addition  P  is  hermitian,  P  -  P‘. 

We  prove  the  identity  of  Eq.  (A  16)  by  considering  the  left-hand  side  of  the  equation.  Now  using 
(AIO) 

( HMH 0~'  -  (HH1  +  HMjH')-'.  (A17) 

Invoking  the  matrix  inversion  lemma  (A2]  we  can  show 

( HH 1  +  HMjH')  '  -  (HHT'  -  (HH'Y'HMjU  +  H'(HH')~]  HMj)~'  H'(HH’)  '  (A18) 

Hence,  using  Eq.  (A9)  it  follows  that 

H'(HH'  +  HMjH')-'H  -  P  -  PMjU  +  PMj)~'P.  (A19) 

Now  consider  the  right-hand  side  of  Eq.  (A16).  Equation  (AI4)  implies  that 

M~'  -  (/  +  PMjP')-'.  (A20) 

Again  using  the  matrix  inversion  lemma,  we  can  show 

(/  +  PMjP‘)-'  -l  -  PMjU  +  P'PMj)-'P'.  (A2I) 

Since  P'P  -  P.  P2  -  P,  and  P  -  P'  then 

PU  +  PMjP')  'P P  -  PMjU  +  PMj)~'P'.  (A22) 

Comparing  Eq.  (A19)  with  Eq.  (A22),  we  see  that  Eq.  (A16)  follows,  and  thus,  (S/A)(—  (5/A)2. 
Hence,  the  configurations  of  Figs.  At  and  A2  are  equivalent  in  that  the  resultant  maximum  signal-to- 
noise  ratios  are  identical. 

This  equivalency  can  be  effectively  utilized  as  follows.  Suppose  there  are  K  external  noise  sources 

with  respective  direction-of-arrival  vectors  (DOAV),  Ak,  and  powers  <r*.  k  —  1 . K.  If  there  is  no 

subarraying  as  seen  in  Fig.  A I  we  can  derive  an  expression  for  the  optimal  signal-to-noise  ratio,  (S/S), 
which  is  a  function  of  Ak,  a\,  k  -  I . K  and  the  steering  vector  S  such  that 

“  F(A\ . Ak,  (t f ,  ...  ,  o-jf,  crfi,  S),  (A23) 

where  F  is  the  functional  relationship  between  the  inputs  and  S/S.  In  order  to  And  the  maximum 
signal-to-noise  ratio  (S/S) |  of  an  array  that  has  a  subarray  transformation  we  use  the  equivalency  of 
the  configuration  in  Fig.  A2  to  show  that 

-  F(PA, . PAk,  <r(* . cr|.  W.  PS).  (A24) 

Equation  (A24)  follows  because  the  (S/S)  as  expressed  by  Eq.  (A23)  holds  for  any  K  arbitrary 
external  input  noise  vectors,  steering  vector,  and  noise  powers.  Hence,  after  the  transformation  by  P, 

the  external  input  noise  vectors  are  PAk ,  k  -  I . K,  the  external  input  noise  powers  are 

unchanged,  and  the  steering  vector  is  PS.  For  example,  if  there  is  only  one  narrowband  noise  source 
with  DOAV,  A\,  and  power  <rj*,  we  can  show  that  if  subarraying  is  nor  used,  then 
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(A25) 


-V  lllsip- 


or 


trj  +  <r?A[At 


IIS'^,111 


where  1 1  - 1 1  denotes  complex  magnitude.  Now  if  we  place  a  subarray  transformation,  H ,  between  the 
antenna  elements  and  the  optimal  weights,  Eq.  (A24)  implies  that  in  order  to  And  the  optimal  subar¬ 
rayed  signal-to-noise  ratio,  we  replace  A\  with  PA\  and  S with  PS  in  Eq.  (A24).  This  results  in 


IlfSlI1- 

o’tf 


W  +  <r\A[PA\ 


IIS'W.II2 


(A26) 


To  derive  (S/N)x  in  Eq.  (A26),  we  used  the  fact  that  PP'  —  P. 
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Appendix  B 

OPTIMAL  2  x  N  BUTLER  SUBARRAY  TRANSFORMATION 


In  this  appendix,  we  find  the  optima)  2  x  \  Butler  subarray  transformation.  The  transformation 
is  optimal  in  the  sense  that  the  signal-to-noise  ratio  loss  that  is  incurred  by  subarraying  using  a  2  x  N 
partition  of  a  Butler  matrix  is  minimized.  We  also  derive  suboptimal  2  x  A'  weighted  Butler  matrix 
transformations  where  it  is  shown  that  subarraying  loss,  although  not  optimal,  is  small  and  within  easily 
computed  bounds. 


Let  us  define  a  subarray  transformation,  H ,  as  follows: 


where 


and 


H-B„ 

(Bl) 

,[11  1 

Vn  i  <>-*■  •  • 

1  1 

(B2) 

.  ,  N-  1. 

(B3) 

Note  that  Bp  is  a  partition  of  a  Butler  matrix.  We  are  given  that  the  N  length  signal  vector,  S,  is 
defined  as 

S-  (l  1  ...  l)r<r|  -  (B4) 

for  an  uniformly  weighted  array,  which  indicates  that  the  mainbeam  is  pointing  in  the  direction  of 
boresight.  In  Eq.  (B4),  o-j  is  the  received  power  of  the  desired  signal  at  each  antenna  element.  Note 
that  multiplying  the  N array  antenna  inputs  by  //results  in  two  channels.  One  channel  will  be  the  main 
channel  in  that  it  is  matched  to  the  steering  vector.  In  the  main  channel,  energy  received  at  an  angle  /9 
off  boresight  will  have  the  following  gain  (or  antenna  pattern): 

g(a)  -  eJko,  (BSa) 

where  a  -  (2nd/\0)  sin/3,  d is  the  linear  array  spacing,  and  A0  is  the  wavelength.  The  auxiliary  chan¬ 
nel  (the  second  channel)  will  have  the  antenna  pattern 

*>>  -  T  (B5b) 

-  g(a  -  *„). 


If  there  is  one  external  jammer  with  power  tr }  and  direction-of-arrival  vector 

A  -  . *<"-'»'♦) r.  (B6) 

where  d  *  (2w<//xo)sin0  and  8  is  the  angle  ofT-boresight  of  the  jammer,  we  would  like  to  know  for 
what  value  of  is  the  output  signal-to-noise  ratio  maximized.  We  will  show  that  #„  should  be  chosen 
as  close  to  the  value  of  d  «s  possible.  Also,  we  show  that  for  arbitrary  N,  choosing  p,  in  this  fashion 


16 


NRL  REPORT  MJ6 


results  in  a  small  signal-to-noise  ratio  loss  from  the  maximum  signal-to-noiae  ratio  that  is  possible  if  no 
subarraying  is  used. 

Let  us  normalize  the  jammer  noise  power  to  the  internal  noise  power,  <r$: 

aj-—-.  (B7) 

o’ 6 

We  use  the  result,  Eq.  (A2S),  of  Appendix  A  to  show  that  the  maximum  signal-to-noise  ratio,  (S/S), 
possible  if  subarraying  is  not  used  is 


5  IISIP 

N  “  (T$ 


S\  Narl 


Oj 

I  IS'A  II2) 

1  +  IMII2<rj 

llsll2  ) 

-  g(<f>)as  we  find  that 

II* 

><*)ll2 

1  +  Ha } 

S 

We  see  that  for  Naj  »  1  and  ||#(4)ll2  «  .  S2  which  is  true  if  the  jammer  is  in  the  sidelobes,  then 

(1)-^ 

If  subarraying  is  used,  then  using  Eq.  (A26),  the  maximum  attainable  signal-to-noise  ratio,  (S/S)wb  is 


"V.  ;  * 

y  ♦ 


5  _  \\PS\\2 

N  sub  0-0 


Qj  IKwi'iP/Oll2 

I  +  \\PA  ||2<x2  \\PS\\ 2 


where 


It  is  straightforward  to  show  that 


P-  H'(HH')~'H  — 


PS  -  S. 

IIWlI2-  Ao-i, 


PA  -  -£*<*>7  +  ^*<#-*)AJ. 


where  1  is  an  A  length  vector  such  that 


1  -  (1  1  ...  1)’ 


and  Af  is  an  S  x  S diagonal  matrix  such  that  A*  -  (A**)  and 

A„  -  e"*-*. 
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Using  Eq.  (BIS)  we  can  show  that 

llP/tll2-  -jp  t||g($)||2A  +  g*(4)g(4  -  +)T'A*T 

+  *(*)*•(* -*)T'a;T 

+  IU(0-^)II2T'a;  aJ] 

•  -^•(llg(^)ll2  +  llg(^  -  *)!l2) 

+  ~2  Ref g*(4)g(4  ~  tl'Mtfr)). 


and  that 


||(PS)'(/><4)||2 


ll^)ll2  +  -^-||^)H2llg(^-*)||2 

+  ~  Re{g*(^)g(i|»)g(^  -  0)}|  •  <rj. 


(B18) 


(B19) 


If  we  set  4  —  4„  where  4„  is  defined  by  Eq.  (B3),  then  g(i>„)  **  0  and  we  can  write  (S/A)Jub  using 
Eqs.  (B14),  (BI8),  and  (B19)  as 


_S )  _  N<rs  . _ llg(*)ll2 _ 

N  sub  W  N&J2  +  llg(^)ll2  +  llg($  -  1 12 


(B20) 


We  see  that  maximizing  (S//V)sub  as  a  function  of  4„  is  equivalent  to  minimizing  the  second  term  of 
Eq.  (B20).  This  term  is  minimized  when  \\g(4  -  d»„>  1 12  is  maximized.  This  occurs  due  to  the  nature 
of  llg(a)||2»  (sin2 Aa/2)/sin2 (a/2)  when  i/i„  is  chosen  to  be  the  nearest  value  of 
n  —  1 . N  —  1,  to  4-  Let  this  value  be  4w 


In  order  to  evaluate  the  signal-to*noise  ratio  loss  incurred  by  subarraying  we  define 


(S/N)Mb 
(S/N)  ’ 


Using  (BIO)  and  (B20)  we  see  that  for  N&}  »  1 

i  _ ik mi _ 


(B2I) 


(B22) 


The  maximum  loss  will  occur  when  4  is  located  equidistant  between  4k  and  #*+|.  This  occurs 
because  llgty  -  *nui()ll2  is  monotonically  decreasing  in  the  region  **_t  <  4  <  where  4k  -  *miK. 
Hence,  since  4  -  4  m,,  can  be  at  most  v/N  and  Liub  decreases  as  llg(4  -  #mM)ll2  decreases,  then  the 
maximum  loss  occurs  when 


Therefore, 


4  4m*» 


w 

N ' 


(B23) 

(B24) 
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In  Eq.  (B24),  we  assumed  that  ir/2N  «  1.  Hence, 


/  max  _  1  _ 

^sub  1 


llg(*)ll2 


llg(^)ll2  +  N&J1 


AN 2 


(B25) 


We  see  that  if  the  jammer  is  in  the  far  sidelobes  then  4N2/ir2  »  ||g(^)||2  and 

/Max  _  j  _  llg(^)ll2 
NcrJ2  + 


AN 2 


(B26) 


Equation  (B26)  implies  that  the  signal-to-noise  ratio  loss  will  be  small.  In  order  to  calculate  the  max¬ 
imum  loss  in  dB  we  use  the  approximation  that  for  small  t 


10  logl0  (!-«)«- 


10 

InlO 


«  dB. 


Hence, 


(/MAX)  10  ikMli  HR 

a“»  )d“  InlO  .  4#>  dB 

N„  +Tr 


10  llg(^)|[2 

ln  10  (Mr2)-1  +  4- 


dB 


where  we  have  normalized  ||g(^)ll2  to  have  a  maximum  value  of  one.  More  explicitly 

1 1£(0)  1 12 : 

For  N<tj  «  A/ it1.  Eq.  (B28)  reduces  to 


j 

jV2  •  2  ^ 
sin2  ir 


<£s«b  x  )<IB  -  10.716  ||g(*)||2dB. 


(B27) 


(B28) 


(B29) 


The  above  results  can  easily  be  extended  to  include  steering  the  mainbeam  off-boresight  and 
weighting  the  antenna  elements  to  obtain  a  given  level  of  sidelobes,  beamwidth,  etc.  However,  the 
adaptive  array  with  no  subarraying  must  also  have  the  same  quiescent  antenna  pattern  determined  by 
the  quiescent  weighting  when  calculating  the  subarray  signal-to-noise  ratio  loss. 

In  addition,  if 

*(*)  -  ak^e**,  (B30) 


where  ak  are  the  quiescent  weights  which  determine  the  antenna  steering,  sidelobe  level,  etc.,  then 
g(+„)  is  not  necessarily  equal  to  zero  and  hence  Eq.(B20)  is  not  exactly  true.  However,  we  can  show 
that  the  terms  that  went  to  zero  in  Eqs.  (B18)  and  (B19)  are  small  compared  to  the  other  terms  so  that 
Eq.  (B20)  is  approximately  equal  to  the  maximum  signal-to-noise  ratio. 

In  any  case,  if  is  chosen  to  be  closest  to  4,  then  Eq.  (B29)  follows  where 

f<*>  -  “  Jjj  «*+.***•  (B31) 
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